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ABSTRACT 

We present broad-band 24, 70 and 160 /im photometry, 5-35 /im and 55-90 /jm spectra of the 
eruptive variable V4332 Sgr from Spitzer observations. The distinguishing feature of the 5-35 /im 
spectrum is an unusually broad absorption feature near 10 /im at the position generally associated 
with silicate-rich dust. Through radiative transfer modeling, we show that this broad feature cannot 
arise from silicates alone but requires the inclusion of alumina (AI2O3) as a dust condensate. The 
case for including AI2O3 is strengthened further by the presence of the AlO radical, a potentially 
important molecule in forming AI2O3. The present detection indicates that porous alumina manifests 
itself through a broadening of the 9.7 /im silicate feature and additionally displays, on the shoulder 
of the silicate feature, a component at ~ 11.5 /im. We discuss how further observations of V4332 Sgr 
may have the potential of verifying some general predictions of the dust condensation process. 
Subject headings: infrared: stars-novae, cataclysmic variables - stars: individual (V4332 Sagittarii) 



1. INTRODUCTION 

Alumina (AI2O3) is considered to play a significant 
role in dust formation around oxygen-rich cool stars. 
Thermodynamic equilibrium calculations indicate that 
it, along with titanium oxides, is one of the earliest 
condensates in the mineralogical cond ensation sequence 
(|Tielenslll990HGail fc Sedlmavrill999D . Observationally, 
there is some debate and uncertainty regarding the spec- 
tral signatures that can be ascribed to alumina that per- 
mit a firm conclusion to be drawn for the presence of 
alumina. In particular features at 11.3 /im and 13 /im, 
seen in the spectra of O-rich AGB and sup ergiant stars 
have often been attributed to alumina (e.g. JSpeck et all 
120001 and references therein). Additionally, various au- 
thors have shown that the inclusion of alumina grains 
in dust models yields better fits to the observed profile 
of the silicate feature at 9.7 /im (especially when the 
feature is broad) and also reproduces better the overall 
infrared spectral energ y distribution (SEP) in selected 
AGB and OH/IR stars dSpeck et alJ l2000b iMaldoni et all 
l2005t IPePew et al.ll2006[ ). In this paper, we present evi- 
dence for alumina dust detection from Spitzer Space Tele- 
scope observations of the nova-like variable V4332 Sgr. 
The distinguishing feature of its mid/far IR spectrum is 
a deep, unusually broad absorption feature at 10 /im. We 
show that this feature cannot be reproduced by silicate 
dust alone and that it is necessary to invoke the presence 
of amorphous alumina grains to explain it. 

V4332 Sgr erupted in 1 994 in what was initi ally consid- 
ered a nova- like outburst ([Martini et al.lll999h . However, 
its subsequent post-outburst evolution to a cool spec- 
tral type indicated that this was not a classical nova 
eruption. The exact nature of V4332 Sgr is of consid- 
erable interest as it, along with V838 Mon and M31 
RV (a red-variable which erupted in M31 in 198 8), may 
form a new class of eru ptive objects (e.g. iMunari et al.l 
l2002tlBond et al.ll2003l ). V4332 Sgr shows an intriguing 
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emission-line spectrum in the optical and near-infrared 
with several rare spectral features. Prominent molecular 
bands o f TiO, ScO, VO and AlO are also seen in the 
optical (jTvlenda et al.l [20051 : iKimesw enger 2006) imply- 
ing an oxygen rich environment. The fundamental band 
of 12 CO at 4.67/im has also been detected in the sourc e 
along with water ice at 3.05/xm (jBaneriee et al.l [20041. 
The IR excess detected in the source, along with the 
molecular and ice features, suggest a cool dusty environ- 
ment around the central star whose effective temperature 
is estimated to be ~ 3250-3280K (jBaneriee et alJl2003t 
iTvlenda et all[2005h . 

2. OBSERVATIONS AND DATA REDUCTION 

V4332 Sgr was imaged with the Multiband Im aging 
Photometer for Spitzer (MIPS: [Rieke et al.ll200l at 24 
and 70 /xm on 15 Oct 2005 and 2 Nov 2006 (70 /im Fine 
and Default modes, respectively). Data at 160 /im were 
obtained on 15 Oct 2005. Spectra were o btained using 
the I nfrared Spectrograph on Spitzer (IRS; iHouck et al.l 
l200l on 18 April 2005 and 19 April 2006. In 2005, low 
resolution (R ~ 60-100) data from ~ 5 — 38 /im and 
high resolution data (R = 600) from ~ 18 — 38 /im were 
obtained. In 2006, high resolution data from ~ 19 — 
38 /im and low resolution data from ~ 5 — 14 /im were 
obtained. In addition, MIPS SED mode data covering 
the wavelength range from ~ 55—90 /im were obtained on 
27 Sept 2005. For the following discussion, data obtained 
in 2005 and 2006 will be referred to as epoch 1 and epoch 
2, respectively. 

The MIP S data were reduced using the Data Analysis 
Tool v3.06 (|Gordon et al.ll2005fh V4332 Sgr was detected 
as a point source by MIPS at 24, 70, and 160 /im, and 
the flux densities were extracted using both PSF fitting 
and aperture photometry. The measured MIPS flux den- 
sities were 2.34±0.07, 1.07±0.11, and 0.12±0.02 Jy at 
24, 70, and 160 /im respectively. At 24 and 70 /im the 
flux densities measured in Epochs 1 and 2 were identi- 
cal within the errors and we report the weighted mean 
of those measurements. The basic instrumental calibra- 
tion of the MIPS SEP mode is similar to that of the 
70 /im imaging mode (|Gordon et al.l [20071 ) but with a 
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Fig. 1. — Top: Epoch 1 (solid curve and points) along with 
the best model fit (dashed curve); see Table HI column 3. Bottom: 
Expanded view around the 10 /im complex. In addition to the best 
fit model, the best fit silicate only model is over plotted (dash- 
dotted curve). It is clearly seen that a pure silicate model yields a 
poor fit to the extended red wing of the data. 

wavelength-dependent illumination correction. The fi- 
nal spectrum of the source is obtained by extracting a 
5-column (49'.'25) aperture from the sky subtracted 2- 
D spectrum after correcting for slit losses. The MIPS 
SED data were scaled by a factor of 0.88 to match the 
flux density in the 70 /im bandpass. Details of the SED 
reduction can be found in Lu et al. (2007, in prep.). 

The IRS data were processed and extracted using the 
Spitzer Science Center (SSC) Pipeline S15.3 product. 
Since V4332 Sgr is a bright point source, the final spec- 
trum of the low resolution modules were combined us- 
ing the SSC post-BCD co-added, nod-subtracted spec- 
tra. No background observation was obtained for the 
high resolution modules and the background level for 
the epoch 1 observations was ~ 0.06 Jy (a factor of 20 
fainter than the source) at ~20 /im and fairly uniform in 
the co-added 2-D low resolution spectra. The SSC post- 
BCD extraction of the high-resolution spectrum agrees 
with the low-resolution spectrum within the uncertain- 
ties. There is no obvious emission/absorption lines seen 
in the high-resolution spectrum; therefore, the final spec- 
tra of V4332 Sgr in epochs 1 and 2 are computed by aver- 
aging both low- and high- resolution modules. In Figure 
[TJ we present the observed IRS and MIPS SED spectra 
of V4332 Sgr and the combined broad band MIPS fluxes 
at 24, 70 and 160 /im. Since there is little apparent 
evolution between the epoch 1 and 2 broad band fluxes, 
we show only the epoch 1 data in Figure [T] Evidence for 
changes in the detailed shape of the SED between epochs 
will be examined below. 

3. RESULTS 

The spectrum of V4332 Sgr is dominated by a deep, 
broad feature at ~ 10 /im, normally associated with the 
presence of amorphous Mg-Fe silicate grains. However, 



this observed 10 /im feature is relatively broad, with an 
additional feature at ~ 11/im and a flattened wing be- 
yond ~ 13 /im. Additionally, signatures of ices and or- 
ganic materials are evident from ^5 — 8 /im (water ice 
at 6 /im, "organics" at 6.8 /im and possibly methane 
ice at 7.7 /zm; see e.g. Bowey & Hofmeister 2005). We 
have modeled the V4332 Sg r spectrum using the radi- 
ation transfer code DUSTY (jlvezic. Nenkova. fc Elitzurl 
11999ft . The limitations of DUSTY include the assump- 
tion of a spherically symmetric shell of material which 
may not be appr opriate for V4332 Sg r as the system 
may have a disk (jBanerjee et al.ll2004ft . However since 
we are interested in exploring the overall shape of the 
observed SED rather than providing a detailed physical 
model of the complete system, we have restricted our- 
selves to the simplest and most generalized assumptions 
in our calculations. As the luminosity of V4332 Sgr is 
poorly known, we have fixed the stellar luminosity for 
V4332 Sgr at 10 4 L Q , the default input value assumed 
by DUSTY. This assumption does not affect the shape 
of the computed spectrum, only the physical scale of the 
system when combined with the dust temperature at the 
inner radius of the shell. We have fit the observed SED 
with two models viz. model 1 that contains silicate dust 
only and model 2 with a mixture of silicate and alu- 
mina dust, where the inclusion of alumina is prompted 
by the presence of a feat ure at ~ 11 urn of t en attributed 
to am orphous alumina (|Sr>eck et al.ll2000l iDePew et al.l 
I2006h . Prompted by the presence of ice absorption at 
~ 6 /im, the grains in both models are coated with an 
ice mantle (20% b y volume). The silica te dust optical 
constants are from iDraine &; Led (119841) while the alu - 
mina optical data used are from IBegemann "et~all (fl997h . 
Corundum was not included in our model as there is no 
evidence in our spectra for the feature at ~13 /zm as- 
socia ted with the presence of this mineral (|DePew et al.l 
2006). We test ed both the ' porous' an d 'compact' alu- 
mina samples of IBegemann et ail (|1997ft ; as there was no 
substantive differences between the models, we restricted 
ourselves to the 'porous' sample for the subsequent mod- 
eling. The ice optical c onstants are those of Wiscombe 
(ftp://climatel.gsfc.nasa.gov/wiscombe/). The range of 
parameters explored is given in Table [TJ The output 
spectra computed using both dust models are shown in 
Figure [JJ It is clearly seen that a pure silicate composi- 
tion matches the observed 10 /im feature poorly. On the 
other hand, the inclusion of alumina in model 2 improves 
the fit significantly. While there is considerable degen- 
eracy in the fits, especially between the optical depth 
and the dust temperature (low temperature, low optical 
depth models are somewhat degenerate with high tem- 
perature, high optical depth models, though they con- 
sistently yield formally worse fits), it is notable that no 
model consisting of only silicate grains provided a satis- 
factory fit to the 10 /im absorption feature. While model 
2 provides a good fit overall fit to the SED and the 10 /im 
feature from 9-12.5 /zm, it reproduces neither the flatten- 
ing beyond ^13 /im, nor the relatively narrow blue wing. 
We explored using diff erent silicate optical constants (e.g. 
lOssenkopf et al.lll992t ) as well as varying the size distri- 
bution (a max and q) but neither approach improved the 
fit in these spectral regions. It is possible that a more 
complex geometry than the simple spherical shell utilized 
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TABLE 1 
DUSTY Modeling 



Parameter 


Value 


Best Fit 


Stellar Luminosity 


10 4 L Q 


fixed 


Stellar Temperature 


3250 K 1 


fixed 


Rout 1 Fiin 


1000 2 


fixed 


Shell p Distribution 


r -2 


fixed 


Composition 


silicates / alumina 


65%/35% 


(fraction by number) 






T9.8 fJLTCL 


2-55 


45 




300 - 1750 K 


1750 K 


Grain Size Distribution 


MRN 3 


fixed 



1 Banerjee et al. 2003. 

2 Maldoni et al. 2005. 

3 Mathis et al. 1977; n(a) ~ a~ q with a m i n max = 0.005, 0.25 fim, 
g = 3.5. 



in DUSTY could improve the fit in these regions. 

The plots in Figure [1] permit a few conclusions to be 
drawn: (i) the substantial improvement in the fits to 
the broad 10 fim feature with the inclusion of alumina 
indicates that alumina is being detected in the source 
and its presence is manifested by a broadening of the 
9. 7 fim silicate featur e. A similar conclusion was reached 
bv lSpeck et all (2000) using data extending to <13.5 /mi. 
(ii) A small, yet clearly discernible, feature is seen at 
11 /im. This feature is attributable to alumina since our 
model calculations show that increasing the percentage 
of alumina in the alumina-silicate mixture of model 2 
enhances the strength of this feature. We note that this 
11/xm feature is seen in a s ignificant number of stars stud- 
ied by ISpeck et alj (|2000r ) implying that alumina grains 
are fairly prevalent. 

4. DISCUSSION 
4.1. The Case for Alumina Condensation 

It is perhaps not surprising to see evidence for alu- 
mina in the dust surrounding V4332 Sgr given the pres- 
ence of the AlO radical in its optical and NIR spectra 
dBaneriee et al.ll2003t iTvlenda et al.l 120051 iKimeswengerl 
2006) and since AlO can play a critical role in the 
production o f alum ina. Laboratory experiments by 
iDemvk et al.l ()2004l ) show that aluminum oxide clusters 
with stoichiomctry A10-(Al203)„ are readily formed in 
laser vaporized metallic Al when quenched in oxygen and 
argon. These clusters are found to be very stable and 
thus ver y good nucleation s i tes fo r dust growth. Addi- 
tionally, |Be!xung^^jbntiin| (|1995h studied the kinetics 
of AlO + O2 reactions at temperatures in the range 300- 
1700K with a view to studying the fate of aluminum in 
combustion. The higher end of this temperature range, 
it may be noted, is very cl ose to the pre dicted condensa- 
tion temperature fl760K;lTielenslll990h of a lumina dust 
around stars. The lBelvung fc Fontiinl (|1995l ) experiment 
shows that AlO becomes oxidized to AIO2 by O2. An ad- 
ditional reaction involving AlO is AlO + O +M = A10 2 
+ M, where the "chaperon" , M, is any atom or molecule 
present that can remove some energy from the newly- 
formed, activated AIO2 (A. Fontijn, private communica- 
tion). Newly formed AIO2 can further interact with AlO 
to generate alumina: AlO + A10 2 + M = A1 2 3 + M. 

These results suggest that AlO is likely to play 
a significant role in the route to AI2O3 formation. 



This conclusion has th eoretical support in the work of 
iGail fc Sedlmavrl (|1999h who show that any possible nu- 
cleation species that can go on to form dust around stars 
should begin with a monomer with exceptionally high 
bond energy. The AlO monomer satisfies this criterion 
and is thus a favored candidate to lead to the forma- 
tion of larger Al m O n clusters that serve as nucleation 
sites for the formation of other grains or to alumina 
grai ns themselves by homo geneous nucleation. While 
the IGail fc Sedlmavrl (|1999f ) analysis is based on ther- 
mal equilibrium considerations, an alternative model is 
the non- equilibrium formati on o f chaotic silicates pro - 
posed by lStencel etall (jl990l ) and lNuth fc Hechtl (|1990h . 
Chaotic silicates form rapidly from a supersaturated va- 
por of metal atoms, SiO, AlO and OH in a hydrogen 
atmosphere (|Stencel et al.l [l990). In the initial stages, 
the higher reduction of Al with respect to Si will lead to 
the preferential formation of Al-0 bonds at the expense 
of Si-0 bonds. This implies that the IR bands of alu- 
mina associated with the Al-0 stretching mode should 
be prominent early in the formation of the chaotic sili- 
cates. However, as the Al atoms become fully oxidized, 
the higher abundance of Si will make the 9.7fim band 
associated with Si-0 bonds dominate. 

Titanium oxides are considered to also be an early 
dust condensate along with alumina. Given that the 
dust that has formed around V4332 Sgr is of fairly re- 
cent origin as implied by the abrupt infrared brighten- 
ing that developed in the source between 2MASS ob- 
servations in 1998 an d subsequent observations in 2003 
(|Baneriee et al.l2003ft . we might expect some signature of 
these species in the sp ectra, though Ti is nearly 30 times 
less abundant than Al (|Speck et a l. 2000) . Bulk titanium 
oxides can have different fo rms: TiO, Ti02, Ti203 and 
T13O5 (|Demvk et all 120041 ) . The most common, Ti0 2 
can exist as brookite and anatase which convert at high 
temperature into rutile which is the third and most com- 
mon form. The rutile spectrum is expected to show a 
broad and strong band at 13-17 fim; the spectrum of 
anatase shows two strong and broad bands around 17 
and 29 fim. The titanium oxide clusters, studied by 
IDemvk et al.l (|2004f ) as possible nucleation sites, have a 
vibrational transition at ~13.5 fim. The above discussion 
pertains mostly to crystalline forms of titanium oxides 
while in V4332 Sgr the amorphous form may more likely 
be present given that both the silicates and alumina are 
of amorphous nature. But a reasonable possibility exists 
that the flattening of the absorption band longward of 
13.5 fim is indicating the presence of titanium oxides. 

4.2. Evolution of the Dust Condensates 

There is potential in the present data to address certain 
aspects of the dust condensation process in astrophysical 
environments. Since the dust formation process in V4332 
Sgr has begun recently - certainly less than 10 years ago 
- and is possibly still ongoing, there are a few spectral 
features that could change with time. As an example, 
in the "chaotic silicates" hypothesis, it is predicted that 
a strengthening of the silicate component of the 9.7 fim 
feature should take place relative to that of the alumina 
and other early condensate components that blend with 
this feature. There is observational support for such evo- 
lution comparing our data between epochs 1 and 2 (see 
Fig[2]). There is a hint that the broad red wing of 9.7 fim 
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Fig. 2. — Enlargement of the 10 fira absorption complex com- 
paring epoch 1 (black) and 2 (blue) data. There is evidence for a 
narrowing of the absorption complex between the two epochs. 

feature has weakened in epoch 2 relative to epoch 1 and 
that there has been an overall narrowing of the 10 //m 
absorption complex. Although the evidence is tentative 
given the small change (~ la) and only two epochs of 
data, such a behavior might be expected as Al and Ti 
atoms become oxidized and Si-0 begins to dominate the 
composition. Further, it is also predicted that the ratio 
of the 10 /im/18 /jm silicate features could be expected 
to change monotonically as silicate dust nucleates and 
anne als in a circumstellar environment (Nu th fc HechtJ 
1990). Freshly nucleated silicates, as laboratory experi- 
ments show, are expected to have a large 10 /um/18 /mi 
ratio i.e. the 18 /jm feature is expected to be weak (con- 
sistent with what is seen in the V4332 Sgr spectrum). 
Thermal processing should increase the strength of the 
18 nm feature. Although the time scales involved in the 
above processes are not clear, it would be worthwhile 
to monitor the spectrum of V4332 Sgr in the future to 
discriminate between possible scenarios for the evolution 
of dust condensates (Dijkstra et al. 2005 and references 
therein) . 

The detection of alumina condensate in V4332 Sgr may 
have implications for the origin of this object. It has 
been postulated that, along with V838 Mon and M31 
RV, V4332 Sgr forms a new class of eruptive variables. 
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The nature of these objects and the source of their erup- 
tive behavior has not been established and ideas ranging 
from an outburst from a compact object in a red giant 
envelope, various thermonuclear events in a single high- 
mass evolved star, a late He-shell flash, stellar mergers 
and even the swallowing of multiple planets by a ma ssive 
star have been proposed fe.g. lTvlenda fc S okcr 2001, and 
references therein). Within the scope of this work, it is 
not possible to discuss in depth the complexities of the 
origin and nature of the V4332 Sgr system. To date, alu- 
mina dust has been almost exclusively detected in AGB 
and other cool evolved stars and, as discussed above, is 
likely a very early condensate in any oxygen rich environ- 
ment, so the detection of alumina in the early condensate 
of V4332 Sgr indicates that conditions in the ejecta are 
similar to those found around cool evolved stars. Thus 
the detection of alumina in V4332 Sgr may provide a con- 
staint on the nature of the eruption if more detailed mod- 
eling of some of the proposed eruption mechanisms rules 
out conditions conducive to the formation of alumina 
grains. In addition, the detection of alumina around 
V4332 Sgr motivates long-t erm monitoring of th e ejecta 
formed around V838 Mon (jBanerjee et al.ll2006h . If in- 
deed these objects are related at a more fundamental 
level than simply having roughly similar outburst char- 
acteristics, we might expect that the conditions in the 
post outburst ejecta of V838 Mon to be similar to those 
in V4332 Sgr. Given that V838 Mon erupted ~8 years 
after V4332 Sgr an d exhibits AlO lines in its spectrum 
(|Lvnch et al.ll2004h . we might detect similar signatures 
of alumina formation around V838 Mon in the coming 
years if both objects do indeed share a common origin. 



Research at the Physical Research Laboratory is 
funded by the Department of Space, Government of In- 
dia. This work is based on observations made with 
the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology under a contract with NASA. Support for this 
work was provided by NASA through Contract Number 
1277253 and 1256424 issued by JPL/Caltech. 



Kimeswenger, S. 2006, Astronomische Nachrichten, 327, 44 
Lynch, D. K., et al. 2004, ApJ, 607, 460 

Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 
Maldoni, M. M., Ireland, T. R., Smith, R. G., & Robinson, G. 2005, 

MNRAS, 362, 872 
Martini, P., Wagner, R. M., Tomaney, A., Rich, R. M., della Valle, 

M., & Hauschildt, P. H. 1999, AJ, 118, 1034 
Munari, U., et al. 2002, A&A, 389, L51 
Nuth, J. A., Ill, & Hecht, J. H. 1990, Ap&SS, 163, 79 
Ossenkopf, V., Henning, T., & Mathis, J. S. 1992, A&A, 261, 567 
Ricke, G. H., et al. 2004, ApJS, 154, 25 

Speck, A. K., Barlow, M. J., Sylvester, R. J., & Hofmeister, A. M. 

2000, A&AS, 146, 437 
Stencel, R. E., Nuth, J. A., Ill, Little-Marenin, I. R., & Little, S. J. 

1990, ApJ, 350, L45 
Tielens, A. G. G. M. 1990, /.From Miras to Planetary Nebulae: 

Which Path for Stellar Evolution?, 186 
Tylenda, R., Crause, L. A., Gorny, S. K., & Schmidt, M. R. 2005, 

A&A, 439, 651 
Tylenda, R., & Soker, N. 2006, A&A, 451, 223 



